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Abstract 

This paper continues earlier work on the quantum evaporation of black holes. This 
work has been concerned with the calculation and understanding of quantum amplitudes 
for final data perturbed slightly away from spherical symmetry on a space-like hypersurface 

at a late Lorentzian time T . For initial data, we take, for simplicity, spherically- 
symmetric asymptotically-flat data for Einstein gravity with a massless scalar field on an 
initial surface E/ at time t = . Together, such boundary data give a quantum analogue 
of classical Einstein/scalar gravitational collapse to a black hole, perhaps starting from a 
diffuse, early-time configuration. Quantum amplitudes are calculated following Feynman's 
approach, by first rotating: T \T\ exp(— i#) into the complex, where < 9 < n/2 , then 
solving the corresponding complex classical boundary-value problem, which is expected 
to be well-posed provided 9 > , and computing its classical Lorentzian action S^iass 
and corresponding semi-classical quantum amplitude, proportional to exp(zS' c i ass ). For a 
locally-supersymmetric Lagrangian, describing supergravity coupled to supermatter, any 
loop corrections will be negligible, provided that the frequencies involved in the boundary 
data are well below the Planck scale. Finally, the Lorentzian amplitude is recovered by 
taking the limit 9 — > 0_|_ of the semi-classical amplitude. In the black-hole case, by studying 
the linearised spin-0 or spin-2 classical solutions in the above (slightly complexified) case, 
for the corresponding classical boundary- value problem with the given perturbative data on 
E p , one can compute an effective energy- momentum tensor < T Ml/ >pff , which has been 
averaged over several wavelengths of the radiation, and which describes the averaged extra 
energy-momentum contribution in the Einstein field equations, due to the perturbations. 
In general, this averaged extra contribution will be spherically symmetric, being of the form 
of a null fluid, describing the radiation (of quantum origin) streaming radially outwards. 
The corresponding space-time metric, in this region containing radially outgoing radiation, 
is of the Vaidya form. This, in turn, justifies the treatment of the adiabatic radial mode 
equations, for spins s = and s = 2 , which is used elsewhere in this work. 

1. Introduction 

This paper is concerned with the problem of finding approximate classical Lorentzian 
(or slightly complexified) solutions of the coupled Einstein gravity/massless-scalar field 
equations, to describe the region of space-time containing outgoing radiation (both spin-0 
and spin-2) in a very large number of modes, generated by quantum-mechanical evapo- 
ration, as a result of nearly-spherical gravitational collapse to a black hole. The space- 
time metric and scalar field <p are split into a 'background' spherically-symmetric 
part (7^ , $), plus perturbations (h^v , 4>^) 1 e tc, which are typically non-spherical. The 
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energy-momentum tensor formed from $ provides the 'matter source' for an exactly spher- 
ical collapse to a black hole. 

In recent papers [1-5], the quantum amplitude for a given perturbative configuration 
(say of the scalar field (p^) on a final hypersurface T,p at a very late time T , was found 
by rotating T slightly into the complex: T — > |T| exp(— i9), for < 9 < 7r/2; then 
calculating the (complex-valued) Lorentzian classical action S c \a, S s for the corresponding 
classical boundary-value problem, which is expected to be well-defined; then computing 
the resulting semi-classical amplitude, proportional to exp(i5' c i ass ); then finally obtaining 
the amplitude for real Lorentzian T by taking the limit of exp(z5 , c i aS s) as 9 — > 0+ . 
Typically, the perturbative scalar-field configuration (p^ given on the late-time surface T, F 
will involve an enormous number of modes, both angular and radial, but with a minute 
coefficient for each such mode. That is, the given <p^> may contain extremely detailed 
angular structure, and also be spread over a considerable radius from the centre of spherical 
symmetry of the background (7^ , $), again with detailed radial structure. Now consider 
the corresponding classical Dirichlet boundary-value problem above, in which one takes 
as given on E F , but chooses </> < ^ 1 - ) | E = (for simplicity), together with the 

complex time-interval-at-infinity T = \T\ exp(— iO) for < 9 < n/2 . The solution for </>W 
will gradually decay towards zero, as one moves from the final surface E^ to earlier times; 
the rate of this exponential decay will be extremely slow when 9 is close to zero. In this 
case, one will find that, at all times t with < t < T (that is, between Ej and E^), the 
classical solution will continue to have complicated angular and radial structure, much as 
does its boundary value ^"^l^ • 

As a result, one must study Lorentzian (or complexified Lorentzian) classical solutions 
for the linearised metric and scalar perturbations (h^} , <p^), which contain classical spin- 
and spin-2 radiation outgoing from the 'gravitational collapse', with detailed structure 
over (typically) an enormous radial extent. The cumulative effective energy-momentum 
tensor, formed quadratically from derivatives of these first-order perturbations, and then 
averaged over several wavelengths of the radiation, so as to produce a smooth averaged 
< T^ u >eff , is expected to be nearly spherically-symmetric, and indeed to have the form 
appropriate to a radially-outgoing null fluid [6]. This viewpoint simplifies enormously the 
description of the 'effective energy-momentum source' due to the wave-like perturbations, 
which then feeds back into the spherically-symmetric background solution (7^ , $), albeit 
over a suitably long time-scale. In this description, the effective energy-momentum contri- 
bution of the emitted radiation can be reduced to just one spherically-symmetric 'density 
of radiation' function of retarded time, instead of an infinite number of multipole or mode 
coefficients for the final boundary data <l>^\-£ • 

The space-metric metric resulting from such a null-fluid effective T^ v is precisely of 
the Vaidya type [7] . This resembles the Schwarzschild geometry, except that the role of the 
Schwarzschild mass M is taken by a mass function m(t,r), which varies extremely slowly 
with respect both to t and to r in the space-time region containing outgoing radiation. In 
this region, the slowly-varying Vaidya metric provides a valid approximation. Of course, 
one does not expect such a relatively simple analytic approximation to the metric and 
scalar field in the strong-field collapse region, where, in the case of a real time-interval T , 
the classical Lorentzian black-hole solution is highly dynamical. 
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In Sec. 2, we discuss the calculation and form of the averaged energy-momentum ten- 
sor < >eff i assuming that both spin-0 perturbations <jft-> and spin-2 (graviton) 
perturbations h^} are present. It is consistently assumed that the time-scale associated 
with typical radiation frequencies is very much less than the time-scale over which the 
background geometry changes. Thus, the wave-like perturbations in the metric and in 
the scalar field can be treated in a WKB approximation, leading to an expression for 
< T^ u >eff ■ One can then verify that this < >eff generates an extremely slow 
evolution of the resulting Vaidya metric. The case of spin-1 (Maxwell) perturbations is 
also discussed. The resulting Vaidya metric in the outgoing-radiation region is described in 
Sec. 3, in different coordinate systems adapted to different aspects of the radiating system; 
this material has also been covered in part in [20] . A brief Conclusion is included in Sec. 4. 

2. High-frequency limit: fields and energy-momentum tensor 

In Sec. 3 of [3], we expanded out the Einstein field equations in powers of e, given 
a perturbative expansion for the classical solution {g^ v , 0) about a spherically-symmetric 
reference or 'background' solution (7^ , $). We write 

g»v{x,e) = llxv {x) +eh$(x) + e 2 h$(x) + ... , (2.1) 

0(x,e) = $(r,r) + e<p {1) (x) + e 2 (p {2) (x) + ... . (2.2) 
At lowest order O(e ), one has the background Einstein and scalar field equations 

<> " \R {0) 1^ = SttTW , (2 .3) 
7^$;^ = , (2.4) 

Here, R$ denotes the Ricci tensor and denotes the Ricci scalar of the background 
geometry 7^ . Covariant differentiation in the background is denoted by a semi-colon 
( ) ;Q or (below) by V a ( ). The background energy- momentum tensor is denoted by 

= - i 7 ^($, a $, /37 -/ 3 ). (2.5) 

The linearised or 0(e 1 ) part of the Einstein field equations reads [9] 



m« a - 2 - 2 n(l)aa - 2 R{o)a ^% 



(2.6) 



where indices on all quantities are raised and lowered using the background metric 7^ . 
As usual [9], we define 

= h$ - \l»„h^ , (2.7) 

where 

= . (2.8) 



3 



Here, R^ua denotes the Riemann tensor of the background geometry 7 Mi , , and T$) de- 
notes the linearisation or 0(e x ) part of the energy-momentum tensor T^ u (x, e), given ex- 
plicitly in Eq. (3.23) of [3]. 

The linearised Einstein equations (2.6) are most easily studied in a 'linearised harmonic 
gauge' [9] in which, by an infinitesimal coordinate transformation, one has arranged that 

hlT = . (2.9) 

At very late Lorentzian times, the background Riemann curvature and the background 
scalar field $ will die off rapidly, whence the linearised Einstein equations (2.6) simplify to 

^ ~ 2 " 2 < } - h {1) ™ + 7m, = , (2.10) 

that is, the linearised vacuum field equations [9], subject also to Eq.(2.9). 

As described in [4], the rate of change with time of the spherically-symmetric back- 
ground geometry 7^ will be extremely small, during the long quasi-static period when the 
rate of emission of quantum radiation by the black hole hardly varies with time. Hence, 
most perturbation modes, for scalar (spin-0) or gravitational (spin-2) oscillations, will be 
'adiabatic' or high-frequency. Within the high-frequency approximation, in a space-time 
without background matter, one may additionally (without loss of generality) impose the 
traceless gauge condition [9]: 

h^ a a = . (2.11) 

In this case, the linearised Einstein field equations (2.6), subject to the transverse-traceless 
(TT) gauge conditions (2.9,11), reduce further [9] to: 

h^ a -2R^ a h^ a = . (2.12) 

At 0(e 2 ), the gravitational field equations give the second-order contribution G^J to 
the Einstein tensor 

G^v = R^u — -Rg^ ; (2-13) 

Following a lengthy calculation, one finds [3] that the Einstein field equations, up to and 
including second order in perturbations, can be written as 

G<2 = SttTW + 8ttT$ + _ G W , (2 .i4) 

Here, T$ is the background energy-momentum tensor (2.5), and G$ denotes — (|)x 

(2) 

the left-hand side of Eq.(2.6). The quantity T^J denotes 

+ \ (h$ - llw h^™ V ff * V p d> 

(2.15) 
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and T' is denned by 



(2.16) 

where the explicit form of C° u will not be needed. 

In the high-frequency limit, after Brill-Hartle or Isaacson averaging [9] over many 
wavelengths (both in space and in time), as summarised in [3], < T^„ > will give the 

leading spin-2 (graviton) contribution to G$ , and Tffi will give the contribution quadratic 
in the scalar fluctuations (p^ . On averaging over many wave periods and over angles, 
following Sec. 3 of [3], one finds for late Lorentzian times in the high-frequency (Isaacson) 
approximation: 

< T' > — — !— < hW - 2 h (1 '>' a h^ a > (217) 

In the transverse-traceless gauge (2.9,11), appropriate for this region of the space-time, 
Eq.(2.3) simplifies to give 

EE < T\ v > TT = ±- < fcW^ h§. v >TT , (2.18) 

The Isaacson averaged energy-momentum contribution of the scalar-field fluctuations, 
taken for simplicity in the late-time region where the background scalar field $ is nearly 
zero, is 

< T$ > = < V^ 1 ) V,</>W > . (2.19) 

Spin-1 Maxwell field perturbations can also be treated in a similar way [16]. For a 
perturbative Maxwell vector potential in the Lorentz gauge 

VMj, 1 ) = , (2.20) 

the Maxwell field equations read [9] 

V^AW - R$ = . (2.21) 
The averaged Maxwell energy-momentum tensor 

<T^>Maxwell = ^- < V? A™ + #° v A^> A^ > Lor (2.22) 



in Lorentz gauge can be simplified using the field equation (2.21) and integration by parts 
in the Isaacson limit, to give 



< >Maxwell = ^ < V M A«« V,A« > Lc 



(2.23) 



For spin-1 Yang-Mills fields, which typically appear when working with locally-supersymmetric| 
theories of supergravity coupled to supermatter [10], a similar but more complicated treat- 
ment can be given. 

Thus Eq.(2.14), averaged over high-frequency fluctuations, and including a spin-1 
Maxwell-field contribution, becomes 



+ i < hW">. tll h^K p , v > TT . 



(2.24) 



Further perturbative corrections to Eq.(2.24) are of a relative size O(e) smaller, and it must 
be understood that one solves Eqs. (2.12,21,24) simultaneously. To ease the notation, we 
henceforth drop the labels Lor and TT. 

For high-frequency (real) massless perturbations (j)( 1 \Aft\h$ , we make an Ansatz 
which is natural for late times: 



A^\x) 



oo £ 



E E 

£=0 m=-£ 

oo £ 

E E E 

£=1 m=-£ P 

oo £ . 

£ £ ^ L 

£=2 m=-£ P JU 



du 



A^ m (t,r,0)e^^' r )/ e + c.c] , (2.25) 
du [(A M )^ mP (t,r,0)e^(*^/ e + cx .] ? ( 2 .26) 

(A^tnP&rW e*^/' + c.c.l ,(2.27) 



where P = + , x denotes the two orthogonal polarisation states for a radially-travelling 
gravitational wave, and also the standard independent polarisation states in the #-and 
(^-directions for a radially-travelling electromagnetic wave [9]. The quantity 6 u {t,r)/e is a 
rapidly-varying real phase, in common to all the spins 0, 1 and 2, which precisely allows 
for the predominantly radial wave- propagation at late times. Defining ip = 0u,(t, r)/e, we 
assume that the first derivative of ip is large in comparison with first derivatives of the 
'amplitude' or of the corresponding tetrad components of (A At ) w £ m p or (A At!y ) w £ m p 

- see below. Schematically, \dip/ip\ 3> |cL4/A| . 

One might expect the high-frequency expansions (2.25-27) for the perturbed scalar, 
Maxwell and graviton fields to give a good approximation to the radially-outgoing radiation 
at late times during the 'long' period when the black hole is radiating in a quasi-static way, 
with its mass 'slowly' decreasing from the initial value Mi . The corresponding approximate 
behaviour of the overall spherically-symmetric 'background' gravitational field at late times 
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is expected to be given by the Vaidya metric [7], as treated in Sec. 3 below. In particular, 
the connection between the late-time high-frequency expansions (2.25-27) and the Vaidya 
metric will be described explicitly in Sec. 3. Of course, in our classical boundary-value 
formulation, with the time-at-infinity T taken to be of the form T = \T\e~ ld , for < 
9 < 7r/2 , the amplitudes A^ m , (A M )^ m p and (A^^e-mP in the late-time high-frequency 
expansions (2.25-27) are related to the scalar, spin-1 Maxwell and spin-2 graviton data on 
the final surface S p , with time T at infinity. Conversely, on following the radiation to the 
past, we reach the strongly-interacting collapse region of the space-time, where both the 
background spherically-symmetric metric 7 Mi , and the scalar field $ may vary rapidly with 
the coordinates t and r . It is only because the background (7^ , $) is the complex solution 
of the boundary- value problem for the spherically-symmetric Einstein/scalar system, with 
a complex time-separation-at-infinity T = |T| exp(— i0), for < 9 < tt/2, that large 
deviations from flatness in the boundary data are expected to be smoothed out in the usual 
elliptic fashion. This is the distinguishing feature of this complex approach; in contrast, 
Lorentzian-signature evolution of the Einstein field equations, including matter, generically 
leads to space-time singularities. If one knew the form of the background solution, then 
(computationally, at least) one could solve the coupled evolution equations for harmonics 
of (say) the perturbed scalar and gravitational fields. One would then, by matching of 
asymptotic expansions [11], have to join the wave-like solutions emerging from the strong- 
field 'collapse region' above onto the high-frequency expansions (2.25-27) for the radiative 
parts of the fields at late times. 

Next, we consider the leading, geometrical-optics limit of the perturbative field equa- 
tions for spin-0, 1, and 2. For spin-0, one has 

7^0^ = , (2.28) 

which is the O^e 1 ) part of the scalar field equation Eq.(2.4), in the late-time limit that 
$ = 0. The spin-1 field equations are given in Eq.(2.21), and the linearised spin-2 field 
equations in Eq.(2.10). Define 

(Mu, = V M W . (2.29) 

Working again in the late-time region and taking $ = there, a straightforward calculation, 
applying the perturbative field equations (2.10,21,28) to the high-frequency expansions 
(2.25-27), together with the TT and Lorentz gauge conditions, leads to the following 
properties: 

, (2.30) 
, (2.31) 
, (2.32) 
, (2.33) 

at lowest order in e . Suppose that the background metric 7^ in this region, written with 
respect to coordinates (t, r, 9, 4>) in the form 

ds 2 = - e 6( *' r ) dt 2 + e a( *' r ) dr 2 + r 2 [d9 2 + sin 2 9 d<£ 2 ) , (2.34) 
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has the Vaidya form [7], but in the coordinate system described in Eq.(3.18) below: 



e-<'.'> = l-^lll , e *(*.D = f-^LVe- . (2.35) 
r \f{m)J 

Here, m(£, r) is a slowly-varying 'mass function', with m = (dm/dt), and where /(m) 
depends on the details of the radiation. Then, for each choice of the integration variable 
uj in Eqs. (2.25-27), Eq.(2.29) has an outgoing- wave solution 

9 u (t,r) = u(t-r*) , (2.36) 

where we define 

f r 



= J dre a ^ , (2.37) 



by analogy with the 'tortoise coordinate' r* = r + 2M ln^(r/2M)— 1 j in the Schwarzschild| 

solution [9,12]. Because of the slowly- varying nature of the background, one has \dtr*\ <C 
1. Note that a general solution of Eq.(2.30): (/c /i ) w (fc A( ) w = 0, would involve a general 
function of (t ± r*). The outgoing-wave solution (2.36) is picked out because we require 
the expansions (2.25-27) to reduce to outgoing Fourier expansions at large radius. 

In a standard fashion, the application of the linearised field equations and gauge 
conditions to the high-frequency expansions (2.25-27) can be carried on to the next order, 
one power of e beyond geometrical optics. For the spin-0 perturbations, one finds 

^V^ + Z^V^ = , (2.38) 

whence 

V M [jA^m| 2 (fc/Jw] = . (2.39) 
For the spin-1 field, one finds 

(A v ) utmP V CT (M^ + 2(k a ) UJ V a (A v ) u£mP = . (2.40) 
Now introduce a polarisation vector (e M )^ m p such that 

(Aft) vim = Ai^emP {e^ujlmP > (2-41) 

(eJutmP (e^Zimp, = Spp' , (2.42) 



A 



lujimP 



{ApUmp {A"y ulmP \ , (2.43) 
where a star denotes complex conjugation. The Lorentz condition implies 

(e^U m p (fc M ) w = . (2.44) 



The Maxwell field equations then imply 



A lu; imp\ 2 (Mc = 0. (2.45) 



Correspondingly, for the spin-2 field (gravitons), one finds 

{A^) ulmP {k a )J a + 2{k a )^{A^) ojlmP , CT = Q . (2.46) 
One then introduces a symmetric polarisation tensor (e^uimP such that 

(A^ u ) u l m p = A2utmP {e^v)ujtmP , (2.47) 

(e^UmP (enUmP> = $PP' , (2-48) 



A 



2ujtmP 



(VW , (2-49) 



where the last equality is valid up to an unimportant phase. The TT condition implies 

= , (2.50) 
7 P (e^)^P = . (2.51) 



A 2 uimp\ 2 {k^u = . (2.52) 



Then the linearised spin-2 field equations imply 
For s = 1 and 2 , write 

AsdraioP = lAsuimup] e s{ *> lmP ; (2.53) 

where o SUJ £ m p is a real phase. From the 'evolution equations' (2.45,52), one finds that 

V„a slvimP = , (2.54) 

provided that 

(fc") w V M (ln|^ mP |) = -Iv"(fc M ) w . (2.55) 

A similar equation holds for the evolution of the spin-0 coefficients A w g m . Now define a 
preferred affine parameter A along the null rays such that 

= , (2-56) 

and the x fJ '(X) are affinely parametrised null geodesies. This can alternatively be written 
in the form 

(M W V(M W = ; (2.57) 
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that is, that the high-frequency waves move along null geodesies. From Eq.(2.55), one sees 
that the amplitude |A sa ,^ m p| decreases if V /x (/c At ) [iJ > , that is, if the null rays diverge. In 
the arguments leading to Eqs. (2.45,53), one finds also that the corresponding polarisation 
tensors are parallely transported along the null geodesic x M (A) . That is, 

V CT (e M )^ mP = , (2.58) 
(k a ) u VaMuimP = . (2.59) 

Of course, such a description in terms of a family of null geodesies will only be valid 
in a comparatively late-time, large-distance region of the space-time. Where space-time 
becomes highly curved, caustics would be expected to develop and the geometrical-optics 
approach would break down. 

Turning again to the Einstein field equations, we calculate the quantities on the right- 
hand side of Eq.(2.24), being the contributions to G^i (7) which are quadratic in the 
spin-0, spin-1, and spin-2 fluctuations. As earlier, <> denotes an Isaacson average over 
times and angles. For an incoherent source of waves, comprising a large number of roughly 
stationary sources (essential for near-spherical symmetry), only terms in Eqs. (2.26,27) with 
£ = £' , m = m' contribute to the average. In the context of Eqs. (2. 26, 27), <> is also an 
average over the random phase 6> w , since a time average. Therefore, at leading order 
0(e~ 2 ) with respect to the high-frequency approximations (2.25-27), one has 

2 f°° 

< V^ 1 ) V^ (1) > = -J2 duJ ( k »)" |A^ m (t,r)| 2 , (2.60) 

e *„ Jo 



em 



and 



2 f°° 

<V lt A a V v A (T > = - du (fc,X (McJ^WmpfcOl . (2.61) 

e tmP Jo 

Further, 

<V (1 / 1 gV^> = - / duJ (M« (Mw \ A 2 ^mp(t,r)\ 2 . (2.62) 

6 JO 



tmP 

Here, we define the quantity \A S0J £mp(t, r)| to be 



2 1 f 2 

\A au i m p(t,r)\ = — / dtt < \A 8U i m p(x)\ >q , (2.63) 

where <><? denotes a time or phase average. Define further 

<T^> = <T$> + <T^> MaX weii +T™ . (2.64) 

Combining the high-frequency approximation with Isaacson averaging leads to the tensor 
(at leading order) 

2 ^ f°° 2 
<T^> = — c s / du (kf^u (k u ) u \A au £ mP (t,r)\ , (2.65) 

e slmP ^° 
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where 

co = 1 , d = -L , c 2 = A- . (2.66) 

As one would expect for a null fluid, one has 

< T a a > =0 (2.67) 

to leading order. The quantity e was regarded as a free parameter above, which helps in 
keeping track of the magnitudes of different quantities in the high-frequency approximation. 
But, given that the terms denoted by exp(z6^(t, r)/e) in Eqs. (2.25-27) are indeed of high 
frequency, one may then set e = 1 in future calculations, without loss of generality. 

One can readily show that < T^ v > transforms as a tensor under (t, r)-dependent 
'background' coordinate transformations. Further, the equations of continuity (2.39,45,52) 
imply the conservation equation V < T^ v >= 0. It is then natural to regard c s |A sw £ m p| 2 ,| 
the total intensity in the high-frequency perturbations, as a measure of the total energy 
density. But while the total energy is independent of the choice of space-like hypersurface, 
the notion of energy density only has significance with respect to a particular space-like 
hypersurface. Denoting the unit future-directed time-like normal vector to the hypersur- 
face by n(°) M , the energy density measured locally by an observer with 4-velocity n^°-* M 
is 

p = n W» < > 

r<x> (2 681 

= 2V / du (n.ku) 2 c s \A 8uem p(t, r)| 2 . 



slmP 



As expected in perturbation theory about the spherically-symmetric background ( i y tiU , $), 
the mass-energy of the massless waves is quadratic in their amplitude, for small \A su ,emp\- 
A further consequence of Eqs. (2.39,45,52) is that the quantity 

N„ = f dSx V^7 (n.k u ) \A soj£mP \ 2 (2.69) 

simP ^ E 

is the conserved total number density (independent of space-like hypersurface) of massless 
waves (massless-scalar particles, photons and gravitons) passing through the space-like 
hypersurface E . 

3. Solution of background field equations 

The Einstein field equations for a spherically-symmetric geometry of Lorentzian sig- 
nature, of the form (2.34), may be derived from the Riemannian field equations, as given 
in Eq.(3.5-ll) of [3], on replacing e b by (-e b ) : 

(1 — e a ) 

a' = 8tt rT rr + ± '- , (3.1) 

r 

V = 8nre a - b T tt - £Llf!2 , (3.2) 

r 

a = 8nrT tr , (3.3) 
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8nT 6 



l-e- a + lre- a (a>-b>) - \ r 2 = &kT 0o = — 2n 

2 2 sin 



(3.4) 



where 

= 



-(l-e- a ) + e~^ a+ V 



d t (ae^ a - b ^ - d r (b', 



\(b-a) 



A 4 

(3.5) 



As usual, for a massless scalar field r), one has T^ v = — \ 7 M i/ ($ )Q ,$ ;i a 7 a/3 ). 

Equations (3.4,5) imply that 



whence, by Eqs.(3.1,2), 



\(a -V) = 



e (a - b) T tt , 
(l-e a ) 



(3.6) 



(3.7) 



We now derive the Vaidya metric, corresponding to a spherically-symmetric null-fluid 
source, in the form (2.34,35) quoted above. Taking the metric form (2.34), we define the 
function m(t, r) by 

e -«(t,r) = 1 _ 2rn{t,r) 



Using Eq.(3.7), we deduce an expression for e b {t, r) : 



1 - 



2m(t,r) 



exp 



dr 



m'(t,f) 



(f-2m(t,f)) 



(3.8) 



(3.9) 



for some f . By elementary flatness at the origin, one must have a — > as r — > on each 
space-like hypersurface. Asymptotic flatness requires setting f = R^ at the outer bound- 
ary, and then taking the limit Roo — > oo . Eq.(3.1) is the Hamiltonian constraint equation 
[9,13]. Use of Eq.(3.8) shows that Eq.(3.1) can be written as a first-order differential 
equation for the mass m(r) inside a radius r at time t = to , say: 



where 



dm 2 

= 47rr P ' 



P = e~ b T tt 



(3.10) 
(3.11) 



is the energy density. Eq.(3.3) is the momentum constraint equation [9,13]. 

We can now determine the background metric at late times, when the energy-momentum| 
tensor is that of the black-hole radiation, following Vaidya [7]. As in Sec. 2, we study the 
gravitational field produced by perturbations whose averaged energy-momentum tensor is 
< >. (Later, we shall move to a coordinate system more suited to retarded radiation.) 
Since the direction (k^)^ of propagation of the radiation in Sec. 2 is null, we choose 



(kr) u e^ b ~^ + (kt) u = 



(3.12) 
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which corresponds to an outgoing- wave boundary condition at large r. Eq.(3.3) also implies 



< T\ > e ^ b ~ a) + <T\> =0 . (3.13) 

The field equations in terms of the metric functions a and b are as in Eqs.(3.1-5), but 

with T^ v replaced by < T^ v >. Using the momentum constraint (3.3), the Hamiltonian 
constraint (3.1), with Eq.(3.13) and p = — < T l t > , we find 

a , + ( e ° ~ !) + he {a-b)/2 = q (314) 

r 

Using Eq.(3.8), one has 

e (b " a)/2 = - = - (kr) u {kt) u e~ a , (3.15) 



that is, 



^ = (^X U-—Y • (3.16) 



Then Eqs.(3.12,15) imply that 



\m' J \ r J 



9 U = e u (m) , (3.17) 



denoting an arbitrary function of m. Finally, one arrives at the Vaidya solution [7,14], in 
the form quoted in Eq. (2.34,35), 

ds 2 = ~ (S)^ 1 -^^)" 1 ^ 2 + (l-^y^ + r 2 dtf , (3.18) 

describing the background space-time 7^ which results from the energy-momentum tensor 
of the high-frequency black-hole radiation. 

A change of variables: (t, r) — > (w, r) can also be found (see below), such that the line 
element [7] is of the Eddington-Finkelstein type [9] : 

ds 2 = _ (1 _ M^A du 2 _ 2dudr + r 2 dti 2 . (3.19) 



Then radially-outgoing null geodesies are precisely paths of constant u . The function m 
is now independent of r and constant along outgoing null rays. In the generic case that 
(dm/du) is not known, it has proved impossible to diagonalise the Vaidya metric and 
to write u as an explicit function of t and r. Since m < and m' = (dm/dr) > 0, 
one finds that, along lines {u = constant}, r increases with increasing t. As in the 
fixed-mass Schwarzschild solution, the Vaidya metric, in the form (3.18), has a coordinate 
singularity where r = 2m(t, r). But, from the (u, r) form (3.19), one can see that the 
apparent singularity in the metric (3.18) at r = 2m(«) is only a coordinate singularity 
[14]. Further [14], the surface {r = 2m(u)} is space-like, lying to the past of the region 
{r > 2m(«)}, since (dm/du) > . In fact, the geometry in the region {r < 2m(u)} (if 
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such a region exists in the 'space-times' considered here, as generated through solution of 
a boundary-value problem) would gradually deviate from the Vaidya form, as one moves 
to the past by (say) reducing u while holding r fixed, since one would reach the region 
of strong-field gravitational collapse. This region can still be described by the diagonal 
metric (2.34), but the full field equations enforce a more complicated coupled solution. 
Provided that the complexified boundary- value problem, outlined in Sec.l, is well-posed 
for a time-separation-at-infinity T = \T\ exp(— i#), for < 9 < n/2 , then the full (complex) 
Einstein/scalar classical solution studied here will be regular at the spatial origin r = 0. 
Indeed, a solution would then be regular everywhere (with respect to suitable coordinate 
charts) in the region between the initial hypersurface Ej and the final hypersurface Hp. 
Since we are considering the case in which both Ej and E^ are diffeomorphic to R 3 , the 
solution should be regular on a region of the form i" x M 3 , where I denotes the closed 
interval [0, |T|]. 

Note further that the regularity of the boundary data (hij , 0), as posed on E/ and 
Ei? , in the spherically-symmetric case, implies that the boundary value m/ j i?(r) obeys 
2m/ ) F(r) < r for all r > . Equality only holds at the centre of symmetry r = . 

We now relate the Vaidya metric, as given in Eq.(3.18), to other coordinate forms of 
the Vaidya geometry. From Eq.(3.7) and from differentiating Eq.(3.15) with respect to r , 
we find 



to" 


to'^ 




2m \ 


2m 






0- 






m' 


TO J 


r J 





(3.20) 



This can be rearranged in the form 



dt(m' (l — ^ 
d r (m> (l- 2 -f)) 



m 
m' 



(3.21) 



which has the solution 

m' (l " ^) = /M , (3-22) 
where /(to) > is arbitrary [7]. Eq.(3.16) can now be rewritten using Eq.(3.22), to give 

e 6(t,r) = e 2*(t,r) ^ _ ^ (3 23) 

where e 2 ^( t,r ^ is defined as 



TO 



2 



e 2^(t,r) = (__) . (3.24) 
\f{m)' 



Hence, the 4-metric can be written in the form 
ds 



2 = - (l- 2m( ^ r) ) dt* + (l- 2 ^ y 1 dr> + r*dtf . (3.25) 
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The Vaidya model in the context of black-hole radiation has, for example, been studied by 
Hiscock [15]. If one chooses f(m) = —m, then in Eq.(3.24) one has tf)(t,r) = 0, whence 
the Vaidya metric of Eq.(3.25) takes a particularly simple form. 

Different choices of the function f(m) correspond to different physical models; once 
f(m) is specified, one determines m as a function of t and r . One can consider the complex- 
ified boundary-value problem in the case (say) that exactly spherically-symmetric initial 
data ('jij , $)/ are given on the initial hypersurface , whereas on the final hypersurface 
Ep , the data consist of a background spherically-symmetric part (7^ , $)f, together with 

weak 'linear-order' fluctuations (hQ, 4>^)f , which correspond in the classical theory to 
emitted gravitons and massless scalar particles. From experience with real elliptic partial 
differential equation theory, one might not unreasonably expect a unique classical solution 
to this Dirichlet boundary-value problem [16,17]. Hence, in particular, in the context of the 
high-frequency approximation of Sec. 2, and of the Vaidya description of the corresponding 
outgoing-wave-source gravitational field of the present Section 3, one would expect that 
(for example) the detailed high-frequency coefficients A w ^ m (t, r, Q), (A M )^ m p(t, r, Q) and 
(^nv)u!£mp(t, r, O) of Eqs. (2.25-27) should be determined by the above Dirichlet boundary 
data. Similarly, the 'free' function f(m) of Eq.(3.22) in the Vaidya description should also 
be determined, and indeed f(m) should be related to the detailed quantities A u e m , etc., 
above. 

At late times and at correspondingly large radii r , the semi-classical mass-loss formula 
should hold to great accuracy: 



where a(m) effectively accounts for the number of particles light enough to be emitted by 
a hole of mass m , and a(m) increases with decreasing m (here allowing for a more general 
model than our Einstein/massless-scalar case). This follows since, as described in Sec. 3 
of [5], the Bogoliubov coefficients for the quantum evaporation of the black hole [17a] are 
given, for practical purposes, by the standard expression 



Here, T SUJ £ m (m) is the transmission probability over the centrifugal barrier of the black hole 
for a mode with spin s, frequency u> and angular quantum numbers (£, m), and m = 2Mu 
is dimensionless, M being the space-like or total ADM (Arnowitt-Deser-Misner) mass of 
the space-time [9]. The original derivation of Eq. (3.27) was in the context where the black- 
hole singularity was taken to persist at late times. But, because of the very-high-frequency 
(adiabatic) method through which the above expression for \(3 su >i m \ was calculated, it 
should still be valid (up to minute corrections) in the case presently being studied, in 
which there is assumed to exist a smooth final boundary Ep with topology M 3 . The 
derivation of Eq.(3.26) then follows as usual. 

In particular, consider the late-time behaviour appropriate to our massless-field model, 
in which a(m) = a = constant. The large-r solution to Eq.(3.26), in the region where 




(3.27) 
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dr 



m > , is of the form 

1. 

m(u) ~ ((M,) 3 + 3a («2-«)) 5 , (3.28) 

where Mj and u% are constants and w. ~ (t — r) at large r . We set m{u2) = Mj for some 
fixed n 2 , so that, as u 2 — > — oo , the space-like (ADM) and null (Bondi) masses [9] agree. 

Now introduce a null coordinate u = u(t, r) , which agrees asymptotically with the 
above requirement u ~ (£ — r) , but which is defined everywhere, via the transformation 

/ 2m(t,r)\ 1 / m , \ m / 2m(t, r)\ 1 

(3.29) 

It may be verified that this definition is integrable, as follows: Using Eq.(3.22), one can 
re- write Eq.(3.29) in the form 

du = " (tr) * - im) dr = ~ (m) = d(M(m)) ■ {3 ' m 

where we define 

It is then straightforward to apply the coordinate transformation implicit in Eq.(3.30), to 
derive the 'null form' (3.19) of the Vaidya metric from the alternative diagonal form (3.18). 

We are now in a position to make further contact with the more detailed treatment 
in Sec. 2 of the high-frequency expansions (2.25-27) for massless spin-0, spin-1, and spin-2 
fields. In the coordinate system (u,r,0,<fi) of Eq.(3.19), we write out Eqs. (2.25-27) in the 
form 

(j>^(u,r,Q) = V / du A^ m ( W ,r,Q)e l ^< r )+c.c. , (3.32) 

etc. The only non-zero component of the null vector (k^)^ is (k r ) LJ [see Eq.(3.34) below], 
which, by Eq.(2.57), is in principle an arbitrary function of u. Thus, the radiation, corre- 
sponding to outgoing waves at null infinity, travels freely along {u = constant} light cones. 
Further, at any point in the Vaidya space-time, a local observer finds only one direction 
in which the radiant energy is flowing. Eq.(2.38) can now be solved to give 

A au £ m p{r,u,U) = , (3.33) 

r 

where h SUJ £ m p(u, CI) is an arbitrary, dimensionless complex function. By this means, the co- 
efficients A saj £ m p can be related to the distribution of weak-field massless-scalar, Maxwell 
and spin-2 graviton data on the final surface E p . 

In the (u, r, 6, (p) coordinate system, the only non-zero component of the Ricci tensor 

is 

Run = - 4 IT = 871 < T ™ > ■ ( 3 - 34 ) 

r z du 
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Then Eq.(2.65) gives 



1 dm 



4nr 2 



-J- = 2 ^2 C * [ duJ [i k u)u;] 2 \A SUJ £mp\ 



simP 

2 



j.2 

slmP 



f 2 

Jo 



2 



(3.35) 



Hence, m{u) can only decrease as u increases; the perturbation amplitudes are non-zero if 
and only if the mass m is changing. Just as the description in Eq.(3.33) of the coefficients 
AsioimP in the high-frequency approximations (2.25-27) leads to a relation between h SLO i m p 
and the perturbative spin-0, 1 and 2 data on the final surface E^, so Eq.(3.35) gives 
(dm/du) and hence m{u) for the 'background part' of the classical solution, in terms of 
the coefficients A S0J e m p • That is, in setting up, as final data for gravitational collapse, 
the 'background part' (7^ , 3>)e f , together with the perturbative part, one should choose 
the radial dependence of the late-time background part 7^ to allow for the mass function 
m{u) corresponding to Eq.(3.33) for the given particle species and spins. 

Provided that the Lorentzian time-interval T at infinity is sufficiently large, one ex- 
pects to study background 3-geometries 7^ on the final hypersurface E p , which are nearly 
flat out to a certain large radius r = R\ , corresponding to the edge of the region in which 
the radiation reaches E p . For r > R± , one expects 7^ to correspond to a slowly- varying 
Vaidya metric, with m(u) = m(T — r) gradually increasing out to a radius R2 which 
corresponds roughly to the beginning of the radiation. At radii r > R2 on the final surface 
Hp , the function m{u) should be approximately equal to Mi , the conserved ADM mass 
of the system. 

At the high-energy end of the emission spectrum, when the black hole approaches 
the Planck scale, the (thermal) mass- loss rate as given by Eq.(3.26) breaks down. The 
amount of energy emitted by the black hole in the final stages of the evaporation will be 
comparable to its mass, to ~ m. To account for the small-mass behaviour of the black 
hole, therefore, the micro-canonical decay rate must be considered. The micro-canonical 
approach is generally more desirable, as the thermal equilibrium between a black hole 
and the exterior radiation is unstable, due to a negative specific heat in the canonical 
ensemble [18]. In addition, there is no information loss in the black-hole evaporation in 
the micro-canonical picture, as energy is conserved. 

For the low-frequency quanta (a; M) characteristic of the majority of the evapora- 
tion process, however, the canonical and micro-canonical ensembles are almost equivalent, 
and one obtains a Planck-like number spectrum and the decay rate Eq.(3.26). The micro- 
canonical decay rate for small m has the form [19] 

— m = 3 f(m) 
^ / Xm p i \ 

~ W p i J (3.36) 
/ m \ 6 / J m 2 \ 

= expl -47T— 2-), 

\m p iJ \ mt) 
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where A is a numerical constant and m v \ , t p i denote the Planck mass and Planck time, 
respectively. This further equation does not have the bad behaviour as m — > 0, associated 
with Eq.(3.26). Indeed, the free function f{m) will naturally have the corresponding form 
for small m, as follows from a dimensional analysis of the field amplitudes. 

4. Conclusion 

We have seen in Sec. 2 how the averaged effective energy- momentum tensor < T^ v > effM 
is calculated, describing, over scales of several radiation wavelengths, the way in which 
wave-like fluctuations in the spin-0 (scalar) field <p^> and the spin-2 (graviton) part of the 

linearised gravitational field h$ contribute quadratically as sources for the 'background' 
spherically-symmetric 4- metric and scalar field (7^,$). A similar description holds for 
the spin-1 Maxwell field. While this contribution is small at any one time, it persists with 
a comparable magnitude during the whole time to during which the black hole radiates. 
Thus, its effects, particularly on the spherically-symmetric background metric 7^, accrue 
secularly; indeed, the averaged contribution < T^ v >eff is precisely such as to determine 
the rate of loss of mass (— rh) in the familiar fashion, leading to the eventual disappear- 
ance of the central concentration of mass, when one works with a complexified time-interval 
T = \T\ exp(— 16), with < 9 < n/2 , for which one expects a classical solution which is 
regular between the initial hypersurface E/ and final hypersurface E p . 

Such an averaged effective energy-momentum source leads to an approximate space- 
time geometry g^ v of the Vaidya type, as described in Sec. 3, valid in the space-time region 
containing the outgoing radiation. This Vaidya description is in turn essential in the 
treatment of adiabatic radial mode equations, as in [3-5]. Subsequently, in [20], we have 
generalised the boundary- value treatment in [1-4], which refers to quantum amplitudes 
with only spin-0 perturbative data on the final hypersurface Tip , at a late time T . In 
[20], we treat the other bosonic cases of spin-1 and spin-2 final data, by means of similar 
but more complicated methods. The fermionic massless spin-i case is treated in [21], 
and a treatment of the remaining fermionic spin-| case is in preparation [22]; this is 
needed as part of the treatment of locally-supersymmetric models. In all these examples, 
understanding of the Vaidya description is essential. 
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